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Abstract 

A detailed calculation of the coherent and incoherent dynamic structure functions of 
the free Fermi gas, starting from their expressions in terms of the one- and semi-diagonal 
two-body density matrices, is derived and discussed. Their behavior and evolution with the 
momentum transfer is analyzed, and particular attention is devoted to the contributions that 
both functions present at negative energies. Finally, an analysis of the energy weighted sum 
rules satisfied by both responses is also performed. Despite of the simplicity of the model, 
some of the conclusions can be extended to realistic systems. 
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Recent measurements of the dynamic structure function S(q,uj) in pure liquid 4 He and 3 He, 
and 3 He- 4 He mixtures at high and low momentum transfer [lj,[^],[[| have revived the interest 
in the theoretical analysis of S(q,u) on both boson and fermion systems , Q ,[§[], It is well 
known that at high momentum transfer q the scattering simplifies and only the incoherent part 
of the response survives. Therefore, splitting S(q,ui) in coherent and incoherent terms may 
help to understand its nature. The main purpose of this letter is to analyze both responses 
for the free Fermi gas at any value of q. In this case they can be calculated exactly because 
the Hamiltonian is easily diagonalized. Despite the simplicity of the model, it can be useful to 
understand certain features of the response in real systems like 3 He. 

The dynamic structure function of a quantum system of ./V identical particles at zero tem- 
perature is defined as the probability of coupling the ground state to all states compatible with 
a momentum and energy transfer (q, uo) after a density fluctuation p q = YljLi e* qrj takes place 

%> ^) = ir y E\(n\p q \0)\ 2 S(E n -E - u) (1) 
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The Fourier transform of S(q,u) generates the density-density correlation factor S(q,t) 
which is the sum of the incoherent and coherent density responses 0. These funct ions are 
defined as 

1 N 

Sinc(q,t) = -J2(e~ l(ir3 e tm e^e- im ) (2) 
i=i 
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S C oh(q,t) = _^( e -^(n-r J ) e -iq- rje im e » q .r Je -im ) (3) 

where H and rj are the Hamiltonian and the position operators, respectively. Notice that with 
the definition given in (||) only those terms where particles i and j are strictly different contribute 
to the coherent response, in contrast to other common definition where S co h(q,uj) is taken as 
the whole density response and thus coincides with ([!]). 

In a realistic interacting system, where the potential entering in H depends on the position 
of the particles, H is non-diagonal neither in momentum nor in configuration space, and the 
evaluation of the expectation values in Eqs. (j|), (||) and (||) becomes a formidable problem. 
Different approximations have been deviced in the past, attempting to discern which are the 



leading contributions in the limits of high ||],|| and low [10] momentum transfer. Even though 
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a considerable amount of information concerning the dynamics of quantum liquids and nuclear 
systems have been gatered using those methods, precise knowledge of the complete spectrum of 
the Hamiltonian is required if one seeks to calculate all three responses exactly. 

One of the few systems for which all eigenvalues and eigenvectors of the Hamiltonian are 
known is the free Fermi gas, where H is diagonal in the momentum representation. In this case, 
the action of the momentum translation operators exp(iq • rj) appearing in Eqs. @ and (||) is 
straightforward 
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and results in an exponential operator whose argument is linear in pj, thus yielding the position 
translation operator of particle j. After projecting on a complete basis of states in configuration 
space, the incoherent and coherent responses become functions of the one- and the semi-diagonal 
two-body density matrices 
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where p is the particle number density, v = q/m is the recoiling velocity of particle j and 



q 2 /2m = mv 2 /2 its associated kinetic energy. Comparing to the 1/q series of the response 



derived by Gersch and coworkers in the early 70's, one readily notices that Si nc (q, t) and S co h{q, t) 
are described by the first terms of the incoherent and coherent expansions, respectively. Higher 
order terms in the series are related to integrals of the potential and vanish when describing the 
response of a free system. As a consequence, the incoherent response of the free Fermi gas is 
entirely given by the Impulse Approximation (IA). 

Due to the close relation between the one-body density matrix pi(x) and the momentum 
distribution n(k), Si nc (q,t) can be Fourier transformed to yield the standard representation of 
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the IA 



S inc (q,uj) 



dkn(fc)^ (k + q)2 - — -u 
(2tt) 3 p7 v M 2m 2m . 



um 



k n(k) dk, 



(7) 
(8) 



4ir 2 pq J\y\ 

Y = mijjjq — q/2 being the West scaling variable and v the degeneracy of each single-particle 
state of definite momentum. The momentum distribution of the free Fermi gas is a Heaviside 
step function n{k) = 9{kp — k) that allows for the occupation of states up to the Fermi surface 
only. Therefore, the integral in Eq. (|8|) can be performed and the result, expressed in terms of a 
new set of dimensionless variables q = qjkp and u = co/ep with ep = k F /2m the Fermi energy, 
is given by 
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which defines the dimensionless incoherent response S-inc^qjUJ^. 

Proceeding in a similar way, the coherent response of the free Fermi gas can be brought to 
a form that closely resembles the IA 
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where n(k, q)) is the generalized momentum distribution introduced by Ristig and Clark 
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Both the semi-diagonal two-body density matrix and the generalized momentum distribu 
tion of the non-interacting system are easily derived 
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where £(z) is the statistical correlation function, characteristic of fermion systems 
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and therefore, inserting n(k, q) ( |l4| ) in Eq.(|lO|), the coherent response may be written as 
v f „. „/ (k + q) 2 k 2 
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The <5(q) term in the integral contributes only at q = 0, while the other carries all the information 
at finite values of the momentum transfer and can be evaluated analytically 
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thus defining the dimensionless coherent response S co h(q,Cj). 

Several conclusions on the behavior of the coherent and incoherent responses can be drawn 
from Equations @ and (|i~7|), some of them being also valid for realistic interacting systems. 

The incoherent response Si nc (q,Cv) of the free Fermi gas is positive defined for all energies 
between q 2 — 2q and q 2 + 2q, and vanishes out of this range. At fixed q, Si nc (q, Cj) is a quadratic 
polynomial in Cv with its maximum located at Co = q 2 , thus being symmetric around this point. 
Hence, both Si nc (q,Cu) and its derivatives are continuous. 

At momentum transfer q < 2, S co h(q,Cv) is non-zero and negative in the range Cv G (q 2 — 
2q,2q — q 2 ). It is split in two different parts, one defined at negative energies and the other 
at positive ones. At fixed q, both functions are quadratic polynomials in Co that differ only 
in the sign of the linear coefficient. This peculiarity gives rise to a symmetric and continuous 
S C oh(q~,&) that, at oj = 0, presents both a minimum and a discontinuity in the first derivative. 
At g's greater than twice the Fermi momentum, the coherent response vanishes. 

The total response S(q,Cj), which is the sum of Si nc (q,Co) @ and S co h(q,Cj) (|l7l), is the well 
known Lindhard function jll]],^] that becomes totally incoherent at q > 2. At g's smaller than 
2, both the coherent and incoherent responses contribute and, even though they are of opposite 
sign in the free Fermi gas, the total response remains always positive. This is a general property 
of the dynamic structure function of all systems, as is apparent from Eq. ([l]). For instance, if at 
a certain energy one of the two responses is negative, its absolute value has to be smaller than 
the value of the other one at that point, which should be positive, as this is the only way to 
produce a total S(q,uj) that is either positive or zero. 



This feature is particularly apparent at negative u/s, as the energy conserving delta appearing 
in the definition of S(q,uj) (Eq. (Q)) forces the non-zero contributions to appear at positive 
energies only. The separation of the response in its coherent and incoherent terms breaks this 
constrain and both Si nc (q,u>) and S co h(q,uj) can locate part of their strength at uo < 0. However, 
as the total response is zero in this range, S co h(q,to) has to be equal to —Si nc (q,Lo) at Q < 0. 
This is once again a requirement for the response derived directly from its definition, and so 
does not depend on the kind of system under study. It is easy to check from equations (9) and 
(|j~7| ) that this holds for the free Fermi gas. 

Figures 1,2 and 3 show the total response of the free Fermi gas and its incoherent and coherent 
parts at q = 0.01, q = 1 and q = 1.9. The upper plots show the Lindhard function (solid line), 
while Si nc (q,(jj) and S co h(q,Co) are drawn below (solid and dashed lines, respectively). There 
are three main regions at q < 2 where the Lindhard function changes its behavior due to the 
different Co dependence of Si nc (q,Co) and S co h(q,Co). At negative Co's lying between q 2 — 2q and 
0, Si nc (q,uj) = —S co h{q,£j) and the total response vanishes. At positive energies smaller than 
2q — q 2 , the addition of the two responses cancel the terms quadratic in Co (see Equations (||) 
and (|T7|)) and the total S(q, Co) is linear. Finally, at q 2 + 2q > Co > 2q — q 2 the coherent response 
vanishes and the dynamic structure function becomes entirely incoherent and quadratic in Co. 
Out of these limits all three functions are zero. 

The different behavior of Si nc (q,Co) and S co h(q,Co) is also reflected in figures 1-3. At very low 
values of the momentum transfer, both Si nc (q, Co) and S co h(q, Co) are large, but strong cancelations 
produce the strength of the total S(q,Co) to be drastically reduced. Notice that in this limit the 
region Co £ (2q — q 2 ,q 2 — 2q) where S(q,Co) has only incoherent contributions covers a range 
ACo = 2q 2 1, which is much smaller than the range where the two responses coexist. As 
a consequence, S(q,Co) looks linear in Co (see figure 1). When q increases, as seen in figure 2, 
both Si nc (q,Co) and S co h{q,Co) are spread and quenched, but the former still shows a maximum 
well centered at Co = q 2 and the latter a minimum at Co = 0. The minimum of the coherent 
response coincides with the point where the two functions defined in Eq. (|l7]) join (Co = 0), 
and the different value between the left and right derivatives at that point produces a peak 
that sharpens with increasing q. At q = 1 the range where the coherent response is non-zero 
becomes maximal, going from Co = — 1 to u) = 1, but the incoherent response has a wider range of 
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existence, and the total S(q,u>) presents two well differentiated regions, one linear and another 
quadratic in G>. As q rises above this value, S co h(q,Gj) reduces both its magnitude and its range 
of definition and, in particular, at q = 1.9 becomes much smaller than Si nc (q,Gi) (see figure 3), 
thus producing a total S(q,Gj) that is mostly incoherent. 

Another useful tool in the study of the dynamic structure function are the energy weighted 
sum rules, which can be also derived for Si nc (q,u) and S co h{q,uj). They are defined to be the 
moments of the different responses 
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and, even though rn^ c (q) and Tnf^ h {q) integrate over all energies, m^ a \q), built from their sum, 
is an equivalent integral of the total S(q, Gi) but obviously restricted to positive energies. 

The first moments read 

"42,(9) = 1 (20) 

mfi(q) = S(q)-1 (21) 

m^M = f (22) 

m^ h {q) = (23) 

and are the same in a correlated system replacing S(q) by the appropriate static structure 
factor. For the free Fermi gas, all odd moments of S co h(q,Gj) cancel due to the symmetry of the 
coherent response around G) = 0, while even ones are related to derivatives of p2- Furthermore, 
as Si nc (q,uj) is entirely given by the impulse approximation, all odd order incoherent moments 
centered at Co = q 2 vanish. Particle correlations in realistic systems may cause the static structure 
factor to be greater or smaller than 1 as a function of q, as happens in liquid helium. Even though 
this change in the sign of rn^ h {q) can be obtained from a S co h(q,uj) symmetric around to = 
that has no nodes, it seems more reasonable to think about a coherent response that has at least 
one point where the function changes sign, located in such a way that both Tn^ h (q) and m^J h (q) 
are fulfilled. Therefore, this difference in the behavior of the zero moment points towards a 
change in the structure of the coherent response entirely induced by particle correlations. 

The contributions of Si nc (q,uj) and S co h(q,uj) at Gj < 0, moments (p0[) to (p3[) differ from 



those obtained integrating over positive energies only. In particular for the Fermi sea, the first 
moments of Si nc {q,uj) integrated at uj > are 

/* OO O P OO / — \ 

^ du> S mc (q, Q) = 1 - - J k(k-~) n(k) dk (24) 

dujujS inc (q,uj) = q 2 + — y k [k - - J n(k)dk (25) 

and coincide with the previous ones at q > 2 where S(q,u>) = Si nc {q,Cj) The residual terms 
appearing in the right hand side of equations (p4j) and ( p5| ) are the manifestation of the negative 
energy contributions, and are the same in realistic systems when the incoherent response is 
described by the impulse approximation. Notice that in such systems the momentum distribution 
n(k) extends up to infinity, and so those terms contribute no matter how large q is. Nevertheless, 
n(k) decreases rapidly with k and they become vanishingly small as q increases. 

In summary, a direct calculation of the coherent and incoherent density responses of the 
free Fermi gas has revealed significant differences between them and with respect to the total 
dynamic structure function. While the latter is known to be non-zero and positive at to > 
only, none of those constrains apply to the coherent and incoherent responses separately. In 
particular, and for the Fermi gas, Scoh(q,uj) is negative and symmetric around uj = and both 
Si nc (q,uj) and S co h{q,uj) present non-negligible contributions at uj < that cancel exactly once 
added up. This cancellation also holds in any realistic system. This behavior is also reflected in 
the energy weighted sum rules. 
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Figure captions 

Figure 1. Comparison between the dimensionless Lindhard Function (solid line above) and the di- 
mensionless incoherent and coherent dynamic structure functions (solid and dashed lines 
below, respectively) at q = 0.01. 

Figure 2. Same as in figure 1 at q = 1. 

Figure 3. Same as in figure 1 at q = 1.9. 
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